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We consider the influence of variation in the fine structure constant a on a promising experiment 
proposed by DeMille et al. to search for variation in the electron-to-proton mass ratio /x using 
diatomic molecules [DeMille et al, Phys. Rev. Lett. 100, 043202 (2008)]. The proposed experiment 
involves spectroscopically probing the splitting between two nearly-degenerate vibrational levels 
supported by different electronic potentials. Here we demonstrate that this splitting may be equally 
or more sensitive to variation in a as to variation in fi. For the anticipated experimental precision, 
this implies that the a variation may not be negligible, as previously assumed, and further suggests 
that the method could serve as a competitive means to search for a variation as well. 
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I. INTRODUCTION 



It is conventionally assumed that the fine structure 
constant a = e 2 /Tic and electron-to-proton mass ratio 
/j = m e /m p are non- varying physical quantities. With a 
successful track record dating back to the developments 
of quantum mechanics and atomic theory, the legitimacy 
of these assumptions is often taken for granted. Modern 
experiments have verified the stability of these respective 
quantities on the fractional level of 10~ 17 and 10~ 14 per 
year 0,0). Still the search for variations in a, \i, or other 
fundamental "constants" has been motivated by theoreti- 
cal attempts to unify gravity with the other fundamental 
forces of nature, with some leading models suggesting 
temporal or spatial dependence of these quantities Q. 
Moreover, astrophysical data has been used to give evi- 
dence for nonzero variation of a and \i over cosmological 
time and distance scales (see also null results of 

Refs. 043) i further motivating efforts to detect signals 
of a and fi variation in the laboratory. 

One promising method to search for variation of /j 
in diatomic molecules has been proposed by DeMille et 
al. (l3| . The envisioned experiment probes the splitting 
between two nearly-degenerate vibrational levels, with 
the two levels being supported by different electronic po- 
tentials. The authors show that for a given variation in 
H, the shift in transition frequency may be large on both 
an absolute scale and relative to the splitting itself. Fur- 
ther still, the authors experimentally identified a favor- 
able transition in 133 Cs2, with one level associated with 
the ground X : 5]+ electronic potential and the other with 
the excited a 3 £+ electronic potential. The and 
a 3 E+ potentials have a common dissociation limit, cor- 
responding to two ground-state Cs atoms. Based on their 
Cs2 analysis, the authors argue that fractional variations 
of /i on the level of 10 -17 or smaller could be detected 
with their technique. 

In this paper we consider the influence of a variation 



on this prospective experiment. Neglect of a variation 
is seemingly justified by the fact that the vibrational 
spectrum is independent of a in the nonrclativistic limit. 
Nevertheless, sensitivity arises from relativistic correc- 
tions to the electronic potential, provoking a closer exam- 
ination of their effects on the proposed Cs2 experiment. 
To facilitate in this analysis, we introduce coefficients 
quantifying the energy shift of individual levels in the 
vibrational spectrum for given fractional variations in fi 
and a; we define these sensitivity coefficients according 
to the relation 
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where v labels the particular vibrational level. We em- 
ploy atomic units (preceding arguments implicitly as- 
sumed this choice), though energies and sensitivity co- 
efficients are often expressed numerically in cm -1 . Fur- 
thermore, we choose the dissociation limit to be our zero 
of energy, with the convenience that this is common to 
both the and a 3 S+ states. These specifications 

are required to unambiguously define the sensitivity co- 
efficients of Eq. (QJ . 



II. EXPRESSIONS FOR SENSITIVITY 
COEFFICIENTS 

Expressions for the sensitivity coefficients of Eq. ([1]) 
may be derived within the framework of the WKB ap- 
proximation. For a given (i, a, and energy E, the phase 
is given by 

r (a,E) 



<f>(H,a,E)= I V2MO0 [E-V(a;r)]dr, (2) 

Jn(a,E) 

where M(fx) oc /i -1 is the reduced mass, V(a;r) is the 
electronic potential, which depends on the fine-structure 
constant in addition to the internuclear separation r, and 
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ri(a,E) and r a (a,E) are the classical inner and outer 
turning points for energy E (for brevity, we will subse- 
quently refrain from writing dependence on p, a, and E 
explicitly). Variations in p, a, and E impart a variation 
in <f> given by 



with the partial derivatives being 
1 




(3) 



(4) 



f dV{r) 
2[E-V(r)} V da 



dr. 



Our choice for the zero of energy ensures both V(r) = 
and dV(r)/da = as r — > oo. 

Energy appears in Eq. ((2J as a continuous variable. 
Boundary conditions imposed on the vibrational wave 
function restrict the values of physically allowed energies, 
and this may be accounted for by the WKB quantization 
condition: 

<P V = (V + §) 7T, (5) 

with the energy satisfying this condition corresponding to 
the vibrational energy level E v (here v, the vibrational 
quantum number, is a non-negative integer). Eq. ([5]) im- 
plies 8cj) v = 0, which combined with Eq. §3§ further gives 



Ap (d<f>/dp) v (d<f>/da) v 
{d<t>JdE) v ^ (d<t>/dE) v ' 



(0) 



where the subscripts v on the right-hand-side denote eval- 
uation at energy E v (evaluation at present-day values of 
a and p is further implied). Taking the definitions 
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the sensitivity coefficients and q a . v appearing in 
Eq. ([T]) may then be associated with g M and q a evalu- 
ated at energy E v . We choose to use the term sensitiv- 
ity factors to distinguish and q a from the sensitivity 
coefficients q^ tV and q a , v - It is worth noting that, al- 
though E v depends on the reduced mass and therefore 
the electron-to-proton mass ratio, the sensitivity factors 
themselves do not. We further note that q M is necessarily 
non-negative. 

The WKB quantization condition ([5]) tells us that as 
we climb the energy spectrum, each new vibrational state 
is associated with an incremental change of n in the 
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FIG. 1. (color online) Potential energy curves for theZ 1 E+ 
and a 3 S+ states of Cs2 as given by analytical functions of 
Ref. [3 and Ref. [l5| , respectively. All units are atomic units. 



phase 4>- We may thus associate p = Tr~ 1 (d<fi/dE) with 
the density of states at energy E. It then follows from 
Eqs. (|4I5I6[) that the sensitivity coefficient q^ v may be 
simply expressed as 



(^+1) , 
2p v 



(9) 



where p v is the density of states at energy E v . 

Expressio n (t9l) was presented in the Letter of De- 
Mille et al. [13]. The authors noted that while the den- 
sity of states is essentially constant for the lowest v (the 
vibrational states being well-described by those of a har- 
monic oscillator), it rapidly increases for the highest v. 
Together with the numerator of Eq. ©, this suggests 
that maximum sensitivity to p, variation occurs within 
the intermediate part of the vibrational spectrum. This 
was a foundational principle for their proposal. While 
Eq. © perhaps provides a more tangible means to visu- 
alize this behavior, Eq. ([7]) will be of greater operational 
use for us. 



III. MODELING V(r) AND dV(r)/da 

Coxon and Hajigeorgiou [l4| and Xie et al. have 
presented analytical potential energy curves based on ac- 
curate fits to experimental data of the X 1 E+ and a 3 E+ 
states of Cs 2 , respectively. These are illustrated in Fig.[T] 
With these curves, we may determine q^ for both states 
directly via Eq. ([7]). To further determine q a [Eq. ([8])], we 
must also know the change in the potentials V(r) with 
respect to a change in a; this we model from computed 
data. 

We begin by describing our determination of dV(r) / da 
for the ground state. We have previously computed the 
potential energy curve for this state in considerin g a vari- 
ation in ultracold atomic collision experiments |l6||; we 
employ this data for our present purposes as well. This 
data was obtained using the relativistic DIRAC compu- 
tational program [l7| within a coupled cluster singles- 
doubles (CCSD) approximation (see Ref. [l6| for further 
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FIG. 2. (color online) Change in the potential energy curve 
of the Cs2 ground state with respect to fractional change in 
a. The circles correspond to computed data. The dashed 
line represents the leading asymptotic term, —a(dCe/da)/r 6 , 
with dCs/da as estimated in Ref. [3]. The solid line is a fit to 
the calculated data points in the shorter range with a smooth 
transition to the asymptotic curve. The inset magnifies the 
region of the maximum. All units are atomic units. 



details). Our data agrees very well with the analytical 
curve in the vicinity of the equilibrium distance and at 
shorter distances, but fails to produce the appropriate 
assymptotic behavior, —C§/r e , ultimately resulting in a 
well-depth 10% too large (see Fig. 1 in Ref. [3). The 
divergence of our data from the correct asymptotic curve 
is due in large, we suspect, to basis set limitations and 
neglect of higher excitations (triples, quadruples, etc.) in 
the CCSD approximation. 

The potential energy curve was computed with vari- 
ous values of a in the neighborhood of a = 1/137. With 
numerical differentiation with respect to a, we obtain 
dV(r)/da. Our computed data is displayed in Fig. [2] In 
principle, an offset for these data points should be chosen 
to fulfill the criteria dV(r)/da = as r — > oo. However, 
as our computed data fails to produce the correct asymp- 
totic behavior for V(r), it is not expected to produce the 
correct asymptotic behavior for dV(r)/da either. For a 
more appropriate representation, we choose to model the 
asymptotic part by —(dCe/da)/r 6 , using the estimate of 
8Cq j da given in Ref. [16[ . We join this asymptotic curve 
smoothly with a fit of our computed data points in the 
shorter range to obtain dV(r)/da; the resulting curve is 
illustrated in Fig. [2] The fractional change in potential 
depth D with respect to fractional change in a is found 
to be d\n(D)/d\n{a) = 0.25. 

A notable feature of dV{r)/da is its non-monotonic 
behavior. This is expected based on two qualitative fea- 
tures of our computed data. Firstly, the potential is 
found to get deeper for an increase in a; correspond- 
ingly, dV{r)/da is negative in the vicinity of r e . Sec- 
ondly, the C% coefficient becomes smaller for an increase 
in a, a result which may be related to the fact that rel- 
ativistic effects diminish the ground state polarizability 
of atomic Cs [l|| [TH, [l9| . The implication is a positive- 
valued dV(r)/da in the asymptotic region. Inevitably 



8V(r)/da must have a maximum in the intermediate re- 
gion, which indeed appears in our computed data. The 
non-monotonic behavior of dV(r)/da has implications 
for q a which will be discussed below. We also note that 
dV(r)/da is essentially linear in the vicinity of r e and at 
shorter distances. 

The a 3 E+ potential energy curve proves more difficult 
to compute accurately than the ground state, and so we 
are led to model dV{r)/da for this state in a less sophis- 
ticated manner. From a relativistic Fock space coupled- 
cluster calculation, we estimate the fractional change in 
potential depth D with respect to fractional change in 
a to be <91n(D)/<91n(cv) w —0.17. This is comparable 
in magnitude to the value obtained for the ground state, 
but with an opposite sign. The negative sign here implies 
that the potential becomes shallower for an increase in 
a. To model 8V(r)/da, we begin by assuming linear be- 
havior in the vicinity of the equilibrium distance and at 
shorter range, taking a fixed value of 0.17(D/a) at r e . 
We smoothly join this to the same asymptotic curve as 
the X : E+ state, with the reasoning that the Cq coeffi- 
cient is common to both states. Qualitatively, we expect 
this to be an accurate depiction of dV(r)/da. Namely, 
this model predicts dV(r)/da to be positive and mono- 
tonically decreasing with r. 



IV. ESTIMATED SENSITIVITY FACTORS 

With V(r) from Refs. 0[H and dV(r)/da modeled 
from computed data, we have the tools required to evalu- 
ate the sensitivity factors and q a via Eqs. (|7l8p . These 
are displayed for both the X 1 S+ and a 3 E+ states in 
Fig. (3Ja) as a function of energy E. The first thing to 
note is that, generally speaking, q^ and q a are of similar 
magnitude for the respective states, with a maximum ab- 
solute value in the range of (0.17 — 0.28)1? in each case. 
However, while q^ peaks at E « —D/A for both states, 
q a peaks in absolute value at the bottom of the potential 
well, E = -D. 

The absolute shift to a transition frequency is deter- 
mined by the difference in sensitivity coefficients for the 
two levels involved in the transition. DeMille et al. sug- 
gested probing the transition between two nearly degen- 
erate vibrational levels in Cs2, with one level supported 
by the X 1 E+ state and the other supported by the a 3 S+ 
state. The larger scale of q^ for the X 1 ^^ state com- 
pared to the a 3 E+ state ensures a significant difference in 
sensitivity coefficients, provided sufficiently bound levels 
are chosen. The specification of nearly degenerate levels 
is motivated by the fact that, in addition to large absolute 
sensitivity of the transition, experimental considerations 
further call for large sensitivity relative to the splitting 
itself a small splitting translates to a large relative 
sensitivity. 

For the described experiment, absolute sensitivity of 
the splitting ui to variations in a and /z may then be 
gauged by Aq^ and Aq a — the differences between sensi- 
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FIG. 3. (color online) Sensitivity factors g M and q a for the X 1 T.f and a 3 E„ states of Cs2 versus energy E. Each curve extends 
from E = to E = —D, where D is the potential depth of the respective state. Panel (a) displays the full extent of the curves, 
while (b) provides a magnification of the low-|_B| region. The vertical dashed line in (b) identifies the approximate location of 
nearly-degenerate vibrational levels (for 133 Cs2), with one level being supported by the ground state potential and the other 
being supported by the excited state potential [13|. All values are expressed in cm -1 . 



tivity factors of the X X E+ and a 3 E+ states — evaluated 
at the energy of the near degeneracy. From inspec- 
tion of Fig. EJa), we see that |Ag M | and \Aq a \ are both 
maximum at the bottom of the a 3 E+ potential well, 
E = —279 cm -1 (Ag M and Aq a lose meaning for energies 
below this). From the perspective of absolute sensitivity 
to fx variation, a transition in this region would be the 
most favorable to probe. However, here the density of 
states is smallest for the two states, making near degen- 
eracies less likely to occur. Approaching the dissociation 
limit (E = 0), the densities rapidly increase and a suit- 
ably small ui is more likely. DeMille et al. have experi- 
mentally identified one transition in 133 Cs2 suitable for 
their proposed method, the near degeneracy occurring at 
approximately E = —17 .4 cm _1 [l3[. Indeed, these levels 
arc high in the vibrational spectrum, having binding en- 
ergy 0.5% and 6% of the potential depth for the X 1 E+ 
and a 3 E+ states, respectively. The actual value of uj 
here depends on the particular transition selected from 
the hyperfme-rotational substructure, DeMille et al. hav- 
ing suggested two possibilities with u> ~ 0.1 cm -1 [l3j . 

Fig. [3Jb) provides a magnification of the sensitivity 
factors in the low-|i?| region, including the location of 
the near degeneracy identified by DeMille et al. In this 
region, we see that the a 3 S+ state has a rather sizable 
contribution to Aq^ and Aq a . Moreover, as q^ is positive 
for both states, these sensitivity factors cancel to some 
degree in Aq^ (i.e., the levels move in the same direction 
with respect to variation in fj,). On the other hand, the 
q a are seen to have opposite signs for the two states, 
thus having a constructive effect in Aq a (i.e., the levels 
move in the opposite directions with respect to variation 



in a). This may be attributed, in large, to the fact that 
a variation in a causes one potential to get deeper and 
the other to get shallower. 

In Table U we present values for the sensitivity fac- 
tors q^ and q a and differential sensitivity factors Ag^ 
and Aq a evaluated at the energies E = —279 cm -1 and 
E = —17.4 cm -1 , corresponding to the bottom of the 
a 3 S+ potential and the location of the near degeneracy 
identified by DeMille et al. For both energies we see that 
\Aq a \ is only a factor of ~3 less than |Ag M |. This ratio 
is found hold for intermediate energies as well. Thus we 
write the following approximate relation for the shift in 
to with respect to variations in // and a: 

\ /z 3 a / 

Specifically for the near degeneracy identified by DeMille 
et al, \K\ 160 cm" 1 . 

In practice, the transition frequency must be measured 
with respect to some reference (clock) frequency uj c , and 
only variation in the ratio Q = ui/u) c may be extracted. 
(We reiterate here that ui and ui c represent the transition 
and clock frequencies given in atomic units; fl is effec- 
tively the transition frequency in units of the clock fre- 
quency.) Evidently the shift in the clock frequency must 
be further taken into account according to the relation 

SQ 5ui 5ui c 

n uj u) c 

DeMille et al. suggest using an optical atomic clock as a 
reference, in which case ui c is essentially independent of [i. 
a dependence has been considered for species currently 
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TABLE I. Sensitivity factors for the X X E+ and a 3 E+ states 
of 133 Cs 2 evaluated at E = -279 cm" 1 and E = -17 .4 cm" 1 , 
these energies corresponding to the bottom of the a 3 E+ po- 
tential and the approximate location of a near-degeneracy be- 
tween vibrational levels [l3|] , respectively. A indicates the dif- 
ference between X 1 E^" and a 3 Ej state values. The density 
of states p is also provided in each case. 



X^S+ a 3 S+ A 



Evaluated at E = -279 cm" 1 
q„ [cm" 1 ] 827 827 

q a [cm" 1 ] -229 47 -276 

p [1/cm" 1 ] 0.065 0.086 

Evaluated at E = —17.4 cm" 1 
<? M [cm" 1 ] 203 48 155 

q a [cm" 1 ] -36 19 -55 

p [1/cm" 1 ] 0.34 0.39 



state and 100% uncertainty to our value of q a for the 
a 3 E+ state. We subsequently conclude that, at the loca- 
tion of the near degeneracy identified by DcMillc et al. 
Ag Q /Ag M = —0.35 with about 50% uncertainty. This 
accuracy is sufficient to draw important qualitative con- 
clusions in the following section. 

Finally, we may compare our result with a cursory es- 
timate found in the review of Flambaum and Kozlov [22| . 
Here the authors predicted a somewhat weaker influence 
of a variation on the proposed Cs2 experiment, effectively 
finding the ratio Aq a / Aq M to be more than a factor of two 
smaller than our present value (the authors also guessed 
the sign to be opposite). Their rudimentary estimate 
used atomic data in place of unkown molecular data and 
neglected important anharmonic effects of the potential; 
our present value is based on a more refined method. 



V. CONCLUSION 



used as optical standards [2(||2l|. For Hg + , \5u) c /u c \ « 
2\5a/a\, while other clocks are much less sensitive to a 
variation. Thus for |-K"|/w 3> 1, as with the transition of 
interest in 133 Cs2, the shift in uj c need not be considered. 

Before concluding, wc briefly discuss the accuracy 
of our estimates, focusing on the results for E = 
— 17.4 cm" 1 . For the X 1 ^ potential, the classical in- 
ner and outer turning points are found to be 6.7 a.u. and 
22 a.u. using the analytical potential of Ref. [T3 | . From 
Fig. [21 we see that dV(r)/da has both negative and posi- 
tive valued segments over this range. We may decompose 
q a into negative and positive contributions accordingly, 
and in doing so we find 

q a = —54 cm" 1 + 18 cm" 1 , 

where the two terms represent the respective signed con- 
tributions. We see that there is a large degree of cancel- 
lation between these contributions, resulting in the value 
q a = —36 cm" 1 found in Table |TJ We tested the sta- 
bility of our result using various models to match our 
data with the asymptotic form of dV(r) /da (e.g., includ- 
ing order-of-magnitudc estimates for variations in C% and 
Cio coefficients) . We find the negative contribution to be 
highly stable with respect to these different models. The 
positive contribution, on the other hand, is found to be 
quite dependent on the overall offset of our data points, 
this being determined by the particular model (see, e.g., 
Fig. [2]). For the alternative models, the positive contri- 
bution never fluctuated by more than a factor of two. 

For the a 3 S+ state, we predict dV(r)/da to be pos- 
itive and to decrease monotonically with r. Here there 
is no cancellation between contributions of different sign 
as with the X 1 E+ state. This, along with the smaller 
scale set by the potential depth, largely validates our less 
sophisticated modeling of dV(r)/da in this case. 

For a rough error estimate, we ascribe 100% uncer- 
tainty to the positive contribution of our q a for the X 1 E+ 



Here we have considered the effect of a variation on a 
prospective experiment to search for variation ofu using 
nearly degenerate vibrational levels in 133 Cs2 [131. We 
estimate this experiment to be only a factor of three less 
sensitive to a variation as it is to p variation. In Ref. [l3| , 
DcMillc et al. argued that this experiment could plausi- 
bly detect variation of p at a fractional level of 10" 17 or 
less. Our result shows that attaining experimental pre- 
cision sensitive to fractional variation of p, at 1 X 10" 17 , 
for example, implies an accompanying sensitivity to frac- 
tional variation of a at 3 x 10 -17 . The most stringent 
laboratory limits to-date allow for annual drift of a at 
this level [l|. Therefore, we conclude that a variation 
may not be negligible for the proposed experiment. 

We can extend further on this conclusion by noting 
that Cs2 was presented in Ref. [HI as a candidate sys- 
tem for a more general experimental method. Ultimately 
other diatomic systems may prove more advantageous, 
and theoretical work is currently underway with the goa l 
of determining optimal systems for this method [23j | . 
Other systems — more specifically, select transitions in 
other systems — may also be significantly more or less sen- 
sitive to a variation than the Cs2 case considered here. 
For example, heavier systems will have a higher density 
of states compared to lighter systems with a similar elec- 
tronic potential. As alluded to earlier, higher density of 
states is a favorable feature for this method as it increases 
the likelihood for near-degeneracies to occur between vi- 
brational levels. At the same time, heavier systems are 
also known to have larger relativistic effects and, presum- 
ably, larger sensitivities to a variation. 

As another example, we note the close relationship be- 
tween the proposal of DeMille et al. and a previous pro- 
posal of Flambaum and Kozlov [24| . Flambaum and Ko- 
zlov suggested using nearly degenerate vibrational levels 
associated with different fine structure components of an 
electronic multiplct, specifically focusing on levels near 
the bottom of the spectrum. While allowing for large 
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relative sensitivity, this particular method does not pro- 
mote enhanced absolute sensitivity. Nevertheless, it is 
rather straightforward to show that such an experiment 
is four times more sensitive to a variation than to [i vari- 
ation, following from the simple a 2 and /i 1 / 2 scaling of 
fine structure and vibrational intervals, respectively. Ev- 
idently the general method of DcMillc et al. has potential 
to be more sensitive to a variation than to /z variation. 

Finally, reinterpreting our above results, we suggest 
that the general method proposed by DeMille et al. to 



probe variation in the electron-to-proton mass ratio in 
diatomic molecules may be an equally viable method to 
probe variation in the fine structure constant. 
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